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ABSTRACT 


n  '1 

The  VV  deactivation  rate  of  h2  by  H20  is  calculated  using  a  hard  core 
repulsive  potential.  The  process  Is  found  to  be  dominated  by  near  resonant  VR 
processes  with  large  changes  in  rotational  motion,  and  the  rate  Is  critically 
dependent  on  the  classical  path  chosen. 

f— 

In  part  2,  a  non-perturbatl  ve  model  of  H20  Is  used  to  calculate  the 
effects  of  centrifugal  distortion  on  absorption  line  Intensities.  The  errors 
are  less  than  lu  percent  for  weak  as  well  as  strong  lines.  The  effects  ot 
collision  broadening  by  h2  and  02  are  presented  using  Anderson  theory  together 
with  alpole  transitions  supplied  by  the  multi -mode  model. 
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A. 


U2-H2U  deactivation  rate 


I.  INTRODUCTION 

Water  vapor  plays  an  important  role  in  the  response  ot  the  atmosphere 
to  infrared  radation,  due  to  its  intense  system  of  absorption  bands.  Water 
also  proves  to  be  an  efficient  absorber  of  vibrational,  rotational,  and  trans¬ 
lational  energy,  and  tends  to  dominate  equil ibriation  processes  such  as  that 
resulting  from  the  propagation  of  laser  radiation  through  the  atmosphere. 
Consequently  it  is  important  to  understand  the  kinetics  as  well  as  the  spectro¬ 
scopy  of  water  vapor. 

* 

In  particular,  tne  deactivation  rate  ot  vibrational ly  excited  h2  by 
h2U  is  of  critical  importance  to  kinetic  cooling.*1*  In  this  process,  a  laser 
pulse  causes  a  transient  cooling  of  tne  atmosphere,  leading  to  thermal  lensing 
rather  than  thermal  blooming  of  the  beam.  The  C02  laser  pulse  induces  transi¬ 
tions  in  atmospneric  C02  from  the  (lou)  to  (uul)  states.  The  ^(uul)  level 
is  roughly  in  resonance  with  N2(l),  and  consequently  most  of  the  C02(UUl) 
energy  is  resonantly  absorbed  by  N2  molecules.  Atmospheric  N2  serves  as  an 
energy  sink  and  the  depleted  C02(luu)  levels  are  filled  by  translational  ener¬ 
gy,  leading  to  a  decrease  in  temperature  or  cooling.  Eventually,  the  N2(l) 
levels  deactivate,  returning  energy  to  translation,  and  heating  overcomes 
cool ing. 
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The  degree  and  duration  ot  the  cooling  process  are  determined  by  the 
deactivation  rate  ot  h^d),  which  is  dominated  by  the  deactivation 

rate.  Thus  it  is  important  to  understand  the  collision  dynamics  of  H^O  and 

we  note  that  experimental  water  vapor  deactivation  rates  are  anomal¬ 
ously  large,  even  when  there  are  no  vibrational  resonances.  In  general,  meas¬ 
ured  values  are  several  orders  ot  magnitude  more  rapid  than  those  predicted  by 
standard  theories  of  VV  exchange.  Since  these  theories  work  well  tor  most 
other  molecules,  it  would  be  interesting  to  know  why  they  fail  for  water  va¬ 
por.  we  nave  calculated  the  collision  rate  for  the  process 

h*U  ♦  M,0  ♦  M-0*.  (1.1) 

c  c  *  i  c 

Our  predicted  rate  is  in  good  agreement  with  recently  measured  values.  In 
particular,  at  Sou0  Whitson  and  f^c Neal  found  K  *  .3  x  10°  sec”^  atm'*,^ 
while  we  predict  K  «  I  x  Iub  sec*^  atm 

II.  COllISlONAl  MECHANISM 

uue  to  its  asymmetry  and  large  rotational  energies,  the  roto-vibra- 
tional  structure  of  water  vapor  is  very  complex,  and  in  fact  there  are  transi¬ 
tions  between  VR  levels  which  are  resonant  with  most  other  molecular  infrared 
transitions.  Conventional  theories  either  neglect  rotational  changes,  or 
assume  that  they  satisfy  dipole  selection  rules.  It  is  a  result  of  these 
assumptions  that  deactivation  rates  are  under-estimated  by  several  orders 

b 


ot  magnitude.  There  is  evidence  that  in  collisons,  large  changes  dt  rotation- 

J  (3  4) 

al  motion  occur,  and  therefore  should  be  taken  into  account.  ’ 

for  our  model  of  the  collision  we  have  chosen  an  interaction 

potential  dominated  by  a  large  repulsive  core  term.  Other  possible  choices 
are  long  ranged  electromagnetic  forces,  and  short  ranged  attractive  chemical 
potentials.  Direct  calculation  shows  that  long  ranged  forces  cannot  account 
for  the  large  rates,  in  spite  of  the  longer  interaction  time.*  Chemical 
forces  are  not  expected  to  be  relevant  to  the  N^-H^O  deactivation.  A  compari¬ 
son  ot  D^O  and  deactivation  rates  supports  the  above  choice.  For  colli¬ 
sions  with  CO,  the  ratio  of  D^Q  and  n^0  deactivation  rates  has  been  measured 


KD_0  -  CO  , 

-  7  x  lu  (2.1) 

H20  -  CO  * 

Because  the  long  ranged  interactions  are  highly  resonant,  they  predict  a  much 
smaller  value  for  the  above  ratio.  On  the  other  hand,  chemical  interactions 
would  imply  a  ratio  of  the  two  rates  to  be  close  to  1,  since  H^O  and  D20  are 
cnemically  identical.  Therefore  we  are  led  to  the  short  ranged  repulsive 
interactions  tor  the  deactivation  of  CO  by  H20.  Since  &2  is  similar  in  struc¬ 
ture  to  CO,  we  likewise  choose  repulsive  interactions  tor  the  n2-H20  deacti¬ 
vation. 

*  A  calculation  was  made  using  the  di pol e-quaorupol  e  interaction.  The  result¬ 
ing  rate  was  tour  orders  of  magnitude  less  than  experiment. 
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The  inter-molecular  interaction  potential  is  assumed  to  be  a  sum  ot 


atonvatom  potentials  of  the  form 


»  ■  Vl>  *  W  *  W  *  W  *  W  *  W 


(2.2) 


where  the  are  atom-atom  Horse  potentials 


Vlrl  ■  0  i'-rU.te-rn‘) 


12.3) 


Horse  potential  forms  were  used  because  they  were  more  convenient  in  determin¬ 
ing  the  classical  path.*0'  The  parameters  U  and  a  were  determined  by  fitting 
the  Horse  potential  to  a  tennaro-Jones  potential,  weighted  over  the  repulsive 

region  ot  interaction.  The  r  in  U.2)  are  interatomic  distances,  and  are  re- 

expressed  in  terms  of  the  separation  of  the  molecular  centers  ot  mass  W,  vibra¬ 
tional,  ano  rotational  coordinates. 

To  specify  the  orientations,  we  choose  the  i ntermol ecul ar  axis  as  the 
polar  axis.  Then  the  m^,u  orientation  is  given  by  two  angles  ( •  *2 ^  ’  and  the 
orientation  by  one  angle  (Sj).  The  use  ot  these  angles  rather  than  coordi¬ 
nates  based  on  a  space-fixed  frame  leads  to  a  helicity  representation  tor  the 
angular  momentum.  This  representation  is  convenient  when  coupling  between 
rotational  and  orbital  angular  momentum  is  negl ected , 1 '  '  which  is  generally  a 

good  approximation.  The  transformation  of  the  Hamiltonian  from  a  space  fixed 

representation  to  the  helicity  representation  introduces  terms  which  couple 
orbital  and  rotational  angular  momentum,  giving  rise  to  a  large  number  of 


coupled  equations.  The  coupling  terms  as  well  as  the  orbital  term  are  propor¬ 
tional  to  Assumng  a  range  tor  these  interactions  to  be  on  the  order  ot 

luA°,  the  effective  time  scales  for  orbital  and  rotational  motions  are 


'UKB  '  80  x  lu'14  s 


(2.4) 


w  -  W/J  x  lu'14  5 


For  most  ot  the  states  important  to  VV  exchange,  J  >  b,  and  r^y  >> 

’i/Kb*  resu^Tmg  in  small  orbital-rotational  coupling.  Thus  we  neglect  the 
coupling  terms.  This  approximation  leads  to  an  unperturbed  Hamiltonian  which 
is  diagonal  in  the  space  ot  states  which  consists  or  the  relative,  vibration¬ 
al,  orbital,  and  rotational  motions. 


In  the  nelicity  representation ,  the  potential  (2.2)  is  of  the  form 


H  MO 

V  ■  V(R.  q  ,  qy  ,  .  «2,  *z) 


(2.5) 


«e  retain  only  the  parts  important  to  relative  motion  and  )  *  HjjlMulU) 


exchange: 


V  ,  a  V  ♦  V, 

’eft  u  1 


U.b ) 


u 


X 


h.,  orientation  angle 


*1 

*2  m  H^O  orientation  angles 

DHh  -  W-K 

a  *  .*4  A° 

In  ( 2 . t> )  we  nave  neglected  k-0  interactions  entirely.  This  is  valid 

to  ten  percent  accuracy  because  the  heavy  oxygen  atom  stays  close  to  the  cen¬ 
ter  or  mass  and  samples  much  less  or  the  potential  surface  than  the  lighter  H 
atoms.  Also,  the  attractive  part  of  the  transition  potential  has  been  ne¬ 
glected  since  the  VR  exchange  is  largely  due  to  the  repulsive  interaction  term. 


The  part  or  U.to)  determines  the  relative  motion  of  the  molecules. 
*e  make  the  classical  path  approx imation ,  which  is  valid  when  the  energy 
defect  at  is  small  relative  to  the  kinetic  energy  t^.  In  our  case  there  are 
contributing  VR  exchanges  with  at  anywhere  from  u  to  buu  cm'*.  Hov*ver,  since 
the  transition  probability  given  by  U.¥)  falls  oft  exponentially  as  a  func¬ 
tion  of  at  with  effective  width  on  the  order  of  i:uu  cm'1,  the  transitions  with 


large  energy  exchange  have  small  weight.  Consequently  the  classical  path 
approximation  should  give  a  fair  picture  of  the  collision  process.  Then  the 


interaction  potential  is 
determined  classically  by 


given  as  a  function  of  time  by  Vj(K(t)),  with  R(t) 
Vu.  The  transition  probability  is  then: 


11 


with  the  matrix  element  <i|V^|f>  is  defined  as: 


.M  .  h 

<  i  I  V 1  f  >  «  <*v  ,  |  V  .(R  ( t) )  | 

Vi'i  1  1 


H..U 
v  t J  f T  t 


(2.  lu) 


The  wave  functions  in  ( 2. lu )  are  for  the  m2o  roto-vibrational  motion,  and  the 
vioration  respectively,  kite  that  both  the  vibrational  and  rotational  dyna¬ 
mics  are  treated  quantun  mechanically  for  m^u,  and  the  matrix  elements  were 

( H  ) 

evaluated  using  our  model.  for  only  the  vibrational  motion  was  deter¬ 
mined  quantun  mechanically  and  was  taken  to  De  a  harmonic  oscillator. 
rotation  has  oeen  neglected  in  the  aDove,  because  its  rotational  energies 
are  small  in  comparison  with  all  other  relevant  energies.  Orbital  motion  is 
taken  into  account  by  the  centrifugal  term  in  the  potential  which  determines 
R(  t). 


Me  note  that  (<!.(*)  is  valid  in  the  limit  that  the  matrix  element 
( 2 . lu )  is  small,  or  single  quantum  exchange  is  valid.  Oirect  calculation 
verifies  that  this  is  a  very  good  approximation  in  the  present  case. 

Although  the  matrix  elements  ( 2. lu )  determine  the  relative  probabili¬ 
ty  of  a  specific  transition,  the  magnitude  of  the  rate  will  depend  critically 
on  the  classical  path,  R(t).  This  is  determined  by  together  with  the  cen¬ 
trifugal  term  i( i*l ) /«^ .  An  ambiguity  arises  since  is  dependent  on  molecul¬ 
ar  orientations  through  the  angles  and  l*nerally,  one  simply  aver¬ 

ages  the  relative  potential 
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(2.11) 


vu  IH>  V  V  ®1)  " 

resulting  in  a  potential  dependent  only  on  R  ot  the  form 

<1|V0|1>  -  t  ( i*l  )/R^  (2.12) 

It  the  molecules  were  rotating  rapidly,  they  would  see  such  an  angle  averaged 
potential,  however,  tor  hard  core  collisions  the  relative  time  scales 

are: 

tcuu  '  2xlu"U  s  (2.13) 

M2U  -14 

fRUT  '  U/Jxl0  S 

hy  .14 

TRUT  *  b4uu/Jxlu  s 

we  theretore  see  that  the  molecules  are  not  rapidly  rotating  during  the  colli¬ 
sion  and  that  the  path  determined  Dy  (2.12)  is  inappropri ate. 

Since  the  molecules  rotate  slowly,  we  use  a  steric  factor  averaging 
technique  to  aetermine  R(t).  In  this  approach  le^,0,*^0)  reter  to  a  particular 
initial  orientation  of  n^O  which  remains  constant  during  the  collision.  After 
the  transition  probability  is  evaluated,  the  result  is  averaged  over  all  possi- 
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Die  values  or  (e^0,#^0).  A  similar  procedure  is  followed  tor  the  molecule. 
It  is  interesting  to  note  that  the  steric  factor  method  or  angle  averaging 
gives  rise  to  transition  probabil ities  larger  by  an  order  or  magnitude  than 
that  obtained  by  potential  averaging.  When  transition  probabilities  are  calcu¬ 
lated  as  functions  of  orientation,  the  angles  which  lead  to  deep  interpenetra¬ 
tion  or  the  molecules  tend  to  dominate  the  average.  Conversely,  when  the 
potential  is  angle  averaged  before  the  calculation,  the  larger  potentials 
associated  with  lesser  penetration  dominate  the  average.  As  Indicated,  the 
difference  in  tne  resulting  rates  can  oe  an  order  of  magnitude  or  more. 

,  The  fact  that  the  method  of  angle  averaging  used  can  affect  the  rate 

by  an  order  of  magnitude  or  more  indicates  that  the  reaction  rates  are  criti¬ 
cally  dependent  on  the  classical  path  chosen. 

The  last  step  in  the  calculation  is  to  average  over  all  Initial 
states.  Relative  momentum  is  floltzman  averaged,  and  orbital  motion  is  aver- 

»  ageo  by  treating  t  as  a  continuous  variable  ana  integrating.^  The  initial 

vibrational  and  rotational  states  are  also  averaged  using  a  boltanan  distribu¬ 
tion. 

Ill.  CUhPOTATIljKAl  RLSULTS 

9  The  results  for  several  different  temperatures  are  compared  with 

experiment  in  Table  I.  uur  values  are  roughly  3  times  that  of  the  most  recent¬ 
ly  measured  values  in  colunn  2,  and  show  the  same  temperature  rise.  Uur  re- 

• 
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suits  agree  to  within  a  factor  of  2  with  colunn  3,  but  do  not  show  the  rapid 
rise  in  temperature  indicated  in  col  win  4.  In  every  case  except  for  colunn  4 
near  room  temperature,  the  predicted  rates  agree  to  within  an  order  of  magni¬ 
tude  of  experiment.  We  flno  much  larger  rates  at  lower  temperatures  than 
those  given  in  reference  (11).  In  Table  11,  we  give  some  calculated  values  at 
atmospheric  temperatures. 
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TAbLE  1.  COMPARISON  OF  CALCULATED  AND  EXPERIMENTAL  RATES 
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IV. 


OIbCUSSlUh 


The  predictions  or  kinetic  cooling  in  reference  (I)  used  deac¬ 

tivation  rates  similar  to  colunn  4  in  Taole  I.  we  predict  a  rate  10  to  1!> 
times  larger  than  colunn  4  at  room  temperature.  Witn  such  a  large  rate,  the 
N^-n^u  deactivation  tenas  to  oe  the  dominant  process  for  tne  relaxation  ot 
vibrational  ly  excited  in  the  atmosphere.  For  example,  the  three  processes 
incluaed  in  reference  (1)  as  effective  in  the  relaxation  ot  h^,*  in  the  atmos- 
pnere  are: 


hOW- 


These  are  roughly  equivalent  in  importance  in  the  deactivation  ot  . 
ever,  it  we  use  our  value  tor  tne  rate  (<t.2c),  the  time  scale  tor  the  latter 

process  pec  ones 

„  *  tt.u  *  lo"4  sec  14.3) 

h^u 

Tnus,  tne  h^u  deactivation  is  an  order  ot  magnitude  more  ettective  than  the 
other  two  processes  considered  in  reterence  (1).  inis  dominance  simplities 
tne  analysis  consioeraoly.  It  turns  out  that  both  the  maximum  temperature 
decrease  ana  the  duration  or  the  kinetic  cooling  became  inversely  proportional 
to  tne  deactivation  rate  alone.  Thus  tne  predictions  ot  these  para¬ 

meters  in  ll)  should  be  decreased  by  a  tactor  ot  lu-lb. 

These  parameters  nave  been  measured  in  an  observation  ot  tne  kinetic 
cooling  ettect  by  Sica1*'.  he  touno  that  in  order  to  explain  his  measure¬ 
ments,  tne  rate  ot  deactivation  at  atmospheric  temperatures  hao  to  be 

U  times  greater  than  tnat  used  in  reterence  (I),  in  agreement  with  our  calcu¬ 
lation. 

v.  Conclusion 

The  calculation  ot  collision  rates  involving  change  ot  vibrational 
quantun  numbers  generally  assunes  reasonant  VV  exchange,  with  an  angle  avera¬ 
ged  potential  determining  the  classical  path.  we  have  found  that  the  large 


N^-H^.0  deactivation  rate  is  due  to  rotational  as  well  as  vibrational  exMange. 
lhe  possibility  or  large  rotational  changes  induced  by  hard  core  collisions 
leads  to  a  very  large  number  of  resonant  exchange  reactions,  tending  to  in¬ 
crease  the  reaction  rate  over  that  ot  pure  vibrational  exchange.  In  addition, 

the  result  was  found  to  be  critically  dependent  on  the  classical  path  chosen. 

Specifically,  steric  factor  angle  averaging  gives  results  larger  by  an  order 
or  magnitude  than  potential  averaging,  and  was  necessary  in  predicting  cor¬ 
rectly  the  order  or  magnitude  or  the  transition  probability.  The  resulting 

rate  is  round  to  be  large  enough  to  dominate  the  deactivation  or  vibrational ly 
excited  in  the  atmosphere.  Consequently  both  the  duration  and  degree  ot 
kinetic  cooling  is  inversely  proportional  to  the  deactivation  rate. 

uur  rate  shows  a  temperature  dependence  wnich  agrees  very  well  with 

that  or  references  (2)  and  (lu).  nowever,  for  T  >  euu  our  rise  agrees  with 
tnat  in  reference  (11).  we  give  the  following  approximation  to  the 
deactivation  rate  as  a  function  of  temperature: 


1 

k  »  *.SS7  x  lu'  e  cm'1  atm1  (tor  T  <  7Su)  (S.l) 

_1 

K  ■  2.a7e  X  lu*  e"**’*1  cm‘l  atm"1  ( for  T  >  75u) 
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8.  H20  BAUD  MUOELS 


I.  Introduction 


To  deal  with  problems  involving  atmospheric  absorbtion  or  emission  of 
radiation,  band  models  describing  the  absorbtivity  ot  individual  molecules  are 
needed.  In  particular,  tor  cases  involving  transmission  ot  high  energy  laser 
beams,  models  are  required  tor  a  wide  range  ot  temperatures  and  pressures. 
Line  by  line  models  satisty  these  requirements. 

The  parameters  which  determine  a  line  by  line  model  can  be  either 
measured  experimental ly,  or  calculated  trom  first  principles.  We  use  a  combi¬ 
nation  ot  the  two  approaches.  The  band  model  parameters  are  determined  from 
the  Multimode  Model 1 d’ tor  water  vapor,  which  in  turn  is  based  on  experi¬ 
mental  quantities.  This  semi-empirical  approach  is  expected  to  be  valuable 
for  several  reasons.  In  the  first  place,  pure  ao-initio  calculations  are 
inaccurate  in  predicting  line  positions  and  strengths.  by  using  a  model  con¬ 
sisting  of  experimental  input  into  a  structure  with  a  theoretical  basis,  we 
can  combine  the  accuracy  of  experiments  with  the  simplicity  of  a  theoretical 
model.  A  model  based  on  experiment  is  expected  to  be  valid  only  in  the  range 
of  temperature  and  pressure  at  which  the  experiments  were  performed,  finally, 
a  semi-empirical  model  can  serve  as  a  guide  for  building  more  accurate  models. 

In  order  to  create  a  line  by  line  band  model,  three  major  components 
are  needed: 
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i.  Line  positions 
Line  strengths 
3.  Line  widths 

The  line  positions  are  determined  directly  from  the  Multimode  Model 
for  -Jater  Vapor.  (8*125  As  is  well  known,  the  detailed  molecular  structure  of 
H9J  presents  unusual  difficulties  as  water  displays  extreme  centrifugal  distor¬ 
tion,  viorational  anhamomcity ,  as  well  as  Fermi  resonances  and  other  types 
of  -node-mode  1 nteract 1 ons .  These  features  dominate  the  molecular  characteris¬ 
tics  of  H  .0,  not  only  qualitatively  but  quantitatively,  and  are  reflected  in 
the  '■oto-vibrat lonal  energy  levels  and  wave  functions.  Consequently  they  will 
s-gni f icantly  affect  the  transition  dipole  moments,  level  positions,  and  line 
intensities.  The  multimode  model  takes  these  features  into  account  and  thus 
should  form  the  basis  for  a  useful  band  model. 

The  line  intensities  are  found  by  evaluating  matrix  elements  of  -he 
dipole  operator  between  wave  functions  generated  by  the  multimode  model.  The 
wave  functions  are  of  the  form: 


*vJtM  “ 


£  i-  iJK*> 

m,K 


d-1) 


Here  C^T  is  determined  by  diagonal  i zi ng  the  quantum  mechanical  Hamil- 
tonion  between  the  basis  states  |m>  and  |JKM>,  where  |m>  is  a  product  of  three 
Morse  oscillator  basis  functions  representing  the  stretch,  bend,  and  asymme¬ 
tric  stretch  modes,  and  |JKM>  is  a  symetric  rotor  wave  function. 


Finally,  the  line  widths  at  atmospheric  pressures  are  assumed  to  be 
dominated  by  collision  broadening.  Thus  one  needs  to  know  how  H.,0  interacts 
with  other  molecules.  We  note  that  collisional  interactions  also  cause  a 
small  shift  to  occur  in  the  line  positions.  Both  the  broadening  and  the 
shifts  can  oe  obtained  from  Anderson* s'  * ^ ^  theory  of  collision  broadening 
wnich  w’ll  b*  discussed  below. 

I l .  Theory 

IIA.  Dipole  Operator 

The  dipole  operator  is  a  function  of  rotational  and  vibrational  coor- 
d;nates.  «e  expand  in  a  Taylor  series  about  the  three  vibrational  coordinates 
q,  ,  q->,  q-j  as  f ol  lows : 


u  *  u000  *  “100ql  *  w010  q2  *  u001  q3  *  w200  ql  *  *** 

•  I,  “tjull  4  *3 

I.J.K 

Here  the  u.jK  all  depend  on  rotational  coordinates  in  the  usual  manner.  Note 
that  the  u^qq.  terms  determine  single  quantum  transitions  in  the 

stretch  mode,  bending  node,  and  asymmetric  stretch  modes,  respectively.  Due 
to  anharmonicity  and  mixing  of  modes  this  separation  is  not  clear,  in  that 
there  will  be  snail  contributions  to  bending  mode  transitions  from  all  ujj^ 
with  (I JK )?  (010) .  This  mixing  complicates  the  determinat ion  of  u.j^  From 
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experimental  quantities,  because  uqqq  is  not  exactly  the  measured  electric 
dipole  moment  and  uq^q  is  not  exactly  the  bending  mode  transition  dipole 
moment . 


The  u,j^  were  fit  to  experimental  quantities  in  two  different  ways. 


Both  methods  involve  truncating  the  series  (2-1)  to  limit  the  number  of  para¬ 
meters.  The  first  method  involves  calculating  a  number  of  individual  line 
strengths,  comparing  with  experimental  line  shifts,  and  varying  the  ujJK  in 
order  to  minimize  the  difference  between  the  two.  The  second  method  assumes 
that  VS  coupling  has  a  small  effect  on  the  integrated  band  intensities.  Ne¬ 
glecting  VS  coupling,  one  can  diagonalize  the  rotational  and  vibrational  parts 
of  the  Hamiltonian  separately.  Taking  only  the  vibrational  wave  functions 

I!  ®:  !->  (2-2) 


|V> 


one  arri ves  at  the  integrated  band  intensity  by  using  a  matrix  element  of 
(2-1)  evaluated  between  purely  vibrational  states  (2-2).  In  terns  of  the 
resulting  matrix  element,  the  integrated  band  intensity  is  just 


„  a  2  ,2 

s°  ■ 


(2-31 


where  v  is  the  frequency  at  band  origin  and  u  is  the  dipole  moment  matrix 
element.  The  fitting  procedure  is  now  to  use  (2-3)  to  calculate  band 
strengths,  compare  with  experiment,  and  vary  the  parameters  utjk  to  minimize 
the  error. 
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The  first  method  is  used  in  the  code  LINE,  and  the  second  in  BAND. 
The  parameters  derived  from  line  fitting  give  intensities  which  are  very 
accurate  in  the  region  of  the  spectrum  in  which  the  lines  were  measured,  but 
less  accurate  elsewhere.  The  band  fitting  method,  though  slightly  less 
accurate  due  to  the  neglect  of  VR  coupling,  gives  results  which  are  consistent 
over  the  entire  spectrum. 

IIB.  Matrix  Elements 


Once  the  dipole  operator  is  determined,  one  needs  to  evaluate  the 
matrix  element  between  VR  wave  functions  such  as  (1-1).  for  simplicity  we 
consider  the  matrix  element  of  a  single  term  in  (2-1), 


i  j  k 
UI JK.  ql  q2  q3 


(2-4) 


denote  by  0 _  the  vibrational  matrix  element 

J  mn 


DmnK  "  <m!  ql  q2  q3 '  n  > 


J  *  J 

and  ErK  the  rotational  matrix  element 


EK'K  *  -  (eR0T)  |  J  K  M> 


(2-5) 


(2-6) 
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The  matrix  element  of  u'  between  VR  wave  functions  is  then 


<v'J't'  |u*|  vJt> 


z 


'm.K' 


jIJK 

mn 


rJ'J  rvJt 
tK'K  LnK 


iJK 


(2-7) 


n!i  n  t  K,K 


The  matrices  0  and  E  are  known  from  the  theories  of  Morse 
oscillators  and  symetric  top  rotor  functions ,  and  the  matrices  C 
arise  f~om  the  Multimode  Model.  One  multiplies  the  matrices  to  yield  each 
element,  and  then  sums  over  all  IJK  to  give  the  total  transition  moment 
between  J  v " J "  r " >  and  |vJt>. 

IIC.  Intensities 


The  line  strengths  are  given  by  the  relation 


S 


-he  v 


9.3Nov|u|2 
- Wc - 


where 


No  -  2.692  x  1019  (TSTp/T) 

| u )  *  dipole  moment  matrix  element  in  esu 

g  ■  line  frequency  in  cm’1 
P  ■  occupation  of  lower  state 


-E-vJt/KT 

Q  (2J ♦  I  )e 

Q  - 


(2-8) 
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IID.  Anderson  Theor 


The  impact  approximation  is  made.  This  assumes  that  most  collisions 


are  bi nary,  i .e. , 


Z  «  1/t, 


(2-9) 


where  Z  is  the  collision  frequency  and  is  the  duration  of  one  collision. 

For  Nt^O,  at  atmospheric  temperature  and  pressure. 


Z  ,  N  d2  /  —  *  109  sec'1 

u  m 


,  ,  1  KT  iri12  -1 

1/rcoii  '  ?  >  sr  '  10  sec 


(2-10) 


One  can  see  from  the  above  for  d  *  3  A0,  N  ^  2.7  x  10",  M  *  10.96  amu,  the 
impact  approximat ion  is  valid. 


Secondly,  the  classical  path  approximation  is  made,  where  it  is 
assumed  that  the  relative  motion  .an  be  described  classically.  This  is  valid 
when  the  energy  exchange  or  energy  defect  is  much  smaller  than  the  kinetic 
energy.  For  the  long  ranged  interactions  which  dominate  collisional  line 
broadening. 


a£  «  100  cm" 


(2-11) 
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|  — I - ■ 


E.  >  200  cm”1 

K  x 

and  the  classical  path  approximation  can  be  used. 

Thus  the  relative  separation  is  assumed  to  be  well  described  by  a 
classical  path,  R(t).  One  generally  assumes  an  addition  that  R(t)  is  essenti¬ 
al  ly  the  straight  1  ine 

R(t)  ■  J  b2  ♦  vV  (2-12) 

wnere  b  is  the  impact  parameter  and  v  the  relative  velocity.  The  interaction 
potential  is  then  implicitly  time  dependent,  v  *  v(R(+)). 


Details  on  how  the  line  broadening  and  line  shifts  result  from  the 
interaction  can  be  found  in  numerous  works.  The  line  broadening  arises  from 
the  virtual  transitions  which  occur  during  the  collision,  and  the  shifts  are 
due  to  changes  in  phase.  The  S  matrix  is  given  by 


S 


-A  (w,b,v) 
e 


e-l/2rM* 


(2-13) 


wnere  A  depends  on  b  and  v  through  the  classical  path  (2-12),  and  is  also  a 
function  of  the  energy  defect  w.  The  real  part  of  A  determines  the  virtual 
transition  probability,  while  the  imaginary  part  determines  the  phase  shift, 
r  is  given  by 
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(2-14) 


rn(b*v)  “ 

in 

Vnm  “  <nlVlm> 

One  readily  sees  from  (2-14)  that  r  is  just  the  sum  of  the  Born  approximation 
probabilities  for  all  the  possible  transitions  out  of  the  state  n,  induced  by 
the  collision  interaction  V.  Since  A(w)  is  a  sun  of  holomorphic  functions  of 
w,  whose  real  part  is  l/2r,  the  imaginery  part  of  A  is  given  essentially  by 
the  Hilbert  transform  of  r(w): 


5 


V  (t) 
nm '  ' 


-iw  t 
nm 


dt  | 


(b.v)  . 


hi-] 


r  (w) 

*  nm v  ' 


nm 


w  -  w 

nm 


dw 


(2-15) 


where 


‘  nm' 


“>  ■  $ 


V  (t)  e"1Wtdt’2 
nm' 


Finally,  one  must  average  over  all  impact  parameters,  all  relative 
velocities,  and  the  populations  of  initial  states.  Assuming  that  r<<l,  i  .e. , 
the  transition  probabilities  per  collision  are  small,  the  line  width  and  the 
line  shifts  for  the  state  n  of  H^Q  perturbed  by  are  given  respectively  by 


*rn  "  bdb  ^  vdv  F(v)  rpl  ‘■(b.v)  P ( J £ ) 


(2-16) 
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(2-17) 


£  n  I  ms  3 


vdv  F(v) 


J1J2 


P(J?) 


where: 

Ji  Jo 

rni  '  is  eq.  (2-14)  for  the  case  N^(J^)  ♦  H20(n)  -  H20(n)  ♦ 

J,Jo 

,nl  c  is  eq.  (2-15)  for  the  case  N^Jj)  ♦  H20(n)  .  H20(n)  ♦  N2(J2) 

N  is  number  density  of  N2 

F(v)  is  Maxwel 1 -Bol tzman  distribution  for  relative  velocities 
P(Jt)  is  probability  for  finding  N2  in  rotation  state  J? 

Note  that  the  initial  states  of  and  relative  motion  are  averaged 
over,  while  the  final  states  are  sunned  over. 

Dividing  through  by  the  factor  c  for  (2-17)  and  he  for  (2-18),  and 
using  the  relation 

N  *  N  ^.16 

o  i 

19  -1 

where  N^  is  lochschnitts  number  (v  2.69  x  10  ),  we  get  the  results  in  cm 

atm"* . 
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Ill.  Calculations  and  Results 


II1A.  Dipole  Moment  Operator 


The  dipole  moment  operator  is  assumed  to  be  of  the  form 


r-*  i  j  k 

u  *  L  <3 

ijk 


where  qi  are  the  vibrational  coordinates  defined  with  respect  to  Figure  1 


t 


When  vibrations  occur,  the  lengths  r^  and  and  the  angle  fluctuate 
about  certain  average  values.  Denoting  the  change  from  rest  positions  as  irlt 
SrT.  s*3.  the  vibrational  coordinates  used  here  are  as  follows: 


(3-2a) 


» 


» 


C  3 -2b ) 


q-  *  46 


(3-2c) 


An  approximate  form  for  (3-1)  can  be  obtained  by  assuming  that  the  dipole 
moment  arises  from  a  charge  inbalance  of  -2q  on  the  oxygen  atom,  and  +q  on 
each  hydrogen  atom.  Then  the  b  and  a  components  of  the  dipole  moment  for  the 
rest  position  can  be  easily  derived: 


Now  Ujj  is  just  the  measured  electric  dipole  moment  of  the  H^O  molecule,  u0* 
Thus  q  can  he  determined  since  the  molecular  geometry  is  known.  Using  (3-2) 
and  tne  definition  of  the  dipole  moment,  it  is  now  possible  to  obtain  the 
relation  (3-1)  to  first  order  in  the  q;'s: 


i 
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(3-4a) 


%0  9 
~  X"  7  uo  q3 


M  +  2m  sin  (e/2) 


ql  ‘  7  tdn  7  uo  q2 


( 3 -4  b) 


Note  that  ua  is  antisymmetric  under  special  inversion,  while  is  symmetric, 
"^i s  conclusion  is  true  for  all  orders  of  q.  Therefore,  the  rotational 
selection  rules  for  and  are  different. 

Using  experimental  values  for  ^  and  6,  ( 3 -4 b )  becomes 

-  1.26  ♦  1 . 86q j  -  i.l9q2  (3-5) 

Note  that  the  signs  of  the  u-  .  terms  vary.  If  the  molecule  was  a  rigid  ro- 

'  J  * 

tor,  with  no  mixing  of  modes  due  to  centrifugal  distortion,  the  signs  of  the 
u)4([  would  be  irrelevant,  since  each  would  be  measured  separately.  In 

reality,  there  is  strong  interference  in  many  cases,  so  that  the  signs  are 
important.  for  fundamental  modes  of  vibration,  if  centrifugal  distortion  is 
present  and  the  sign  of  u, ^  is  positive,  the  R  branch  is  strengthened  and  the 
?  branch  is  weakened.  The  reverse  is  true  for  u.  ^  negative.  In  the  case  of 
hot  bands  the  distortion  effects  are  more  complicated. 

Equation  (3-5)  rperesents  a  rough  approximat ion  to  the  dipole  opera¬ 
tor.  As  stated  in  Section  II,  we  have  determined  the  either  by  fitting 

to  individual  line  strengths,  or  to  integrated  band  intensities. 
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In  particular,  we  consider  as  an  example  of  our  model  the  interaction 
ot  the  modes  (uuu),  (ulu),  (U2U),  ( 1 UO ) .  (Ilu),  (12U),  and  ( 00 1 ) .  The  experi¬ 
mental  band  strengths  are  given  below  in  Table  3: 


TABLE  3.  BAND  STRENGTHS  FOR  H,0  (16) 


The  exact  form  of  the  dipole  moment  operator  considered  was 


(3-6 1 


u000  *  w100ql 


w010q2  * 


u001q3  + 


u 1 10q 1 q2  +  u020q2  +  w120qlq2 


This  gave  7  parameters  for  12  pieces  of  data,  and  the  results  of  the  fitting 
procedures  for  this  simple  case  are  given  in  Table  4, 


TABLE  4.  EMPIRICAL  DIPOLE  MOMENT  OPERATOR 


(Debyes) 

u000 

1.851887 

u100 

.457962 

u010 

-.781140 

U1 10 

-.497356 

w  120 

-.789012 

u020 

.014259 

w001 

1.4277 

Note  that  and  uq,q  agree  in  order  „f  magnitude  with  (3-5),  and  the  signs 

are  the  same.  Thus  one  can  understand  the  deviation  of  the  Ho0  absorption 
spectrin  from  a  rigid  rotor  spectrum  qualitatively  in  terms  of  the  simple 
geometric  nodel  of  F^O  shown  in  Figure  1. 
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1 1  IB.  Intensities 


The  intensity  for  a  particular  transition  is  given  by  the  following 

relation: 


1  "  N0P(TSTp/T)(l  -  e 


FT”  3»  r| 


-?  -1 
.2  cm  atm  1 


(3-7) 


where  Nq  ■■  Lochschmidt' s  number 

_  1 

r  *  line  frequency  in  cm 

7 

|  uT  *  dipole  moment 
TSTp  »  273. lo°K 

?  *  Boltzman  distribution  function  for  lower  state 


-  q  UM)e‘£jT/KT  TK 


we  also  give  the  dipole  moment  matrix  element. 


S  -  I u| 2  (2J+1) 

2  2 
| u |  in  units  of  (debyes) 


(3-3) 


Given  S  one  can  calculate  the  line  intensity  for  any  temperature  using  (3-7). 

S  is  determined  by  multiplying  through  the  matrices  C*.  D,  C,  and  E 
as  in  (2-7).  C  represents  the  VR  wave  function,  D  the  matrix  of  the  q's,  and 


L  is  the  matrix  ot  the  rotational  coordinates.  We  note  that  since  i^O  belongs 
to  the  point  symmetry  group  (^y,  ma*r’x  ^  factorizes  Into  4  blocks.  If 
this  is  taken  into  account,  the  multiplication  process  can  be  speeded  up  by  a 
factor  ot  (4 )*■  *  lb. 


In  Tables  b,  6,  and  7,  we  list  the  matrix  elements  S  in  (3-tt)  for 
various  transitions  in  the  (ulu)  band  ot  water  vapor.  The  results  are  com¬ 
pared  witn  experimental  measurements  of  individual  line  strengths,*17'  results 
due  to  Benedict, (l0>  and  the  "rigid  rotor"  line  strengths.  The  latter  are 
calculated  by  neglecting  the  rotational  mixing  ot  modes  due  to  centrifugal 
distortion.  The  “F  factor”  is  obtained  by  dividing  the  real,  centri f ugal ly  dis¬ 
torted  values  by  the  rigid  rotor  values.  The  F  factors,  which  are  related  to 
the  centrifugal  distortion,  are  small  near  the  band  centers,  and  grow  larger 
very  rapidly  as  the  distance  from  the  band  center  increases.  The  original 

F-factor  calculation  assumed  mixing  only  between  the  fundamental  and  ground 
state,  while  in  our  model  mixing  occurs  between  all  the  mooes.  For  regions  of 
the  spectrun  given  in  Table  6,  the  differences  are  not  significant,  but  in  the 
extreme  wings  of  the  spectrum  there  are  order  of  magnitude  differences.  Line 

strengths  for  the  region  of  22uu  cm  1  in  the  (Ulu)  band,  are  given  in  Table 

-3 

7.  Note  that  the  F  factors  there  are  as  small  as  lu  .  Only  a  multimode 
model  can  accurately  predict  the  intensities  in  these  extreme  regions. 
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V 

lower 

r-] 

T 

upper 

V 

exp 

(x  lu'4) 

s 

calc 
(x  lu“4 

SB 

(x  10-4) 

' 

^rlgid 
(x  lu-4) 

7  „ 

0  . 

1455.2b 

899 

873 

1001 

1019 

•0 

-0 

b-2 

4-2 

14t>V.  27 

bub 

441 

517 

411 

i  31 

*.i 

14o4.*3 

21b 

2ul 

267 

15b 

4-l 

J-i 

14/2. ub 

389 

370 

443 

335 

b-b 

5’5 

147 j. 52 

808 

730 

846 

645 

Jo 

2o 

14«7.3b 

3S2 

320 

391 

285 

°-3 

°-5 

UbV.bS 

5ut 

481 

611 

320 

5  , 

4  . 

14*6. 2b 

600 

5b7 

668 

669 

-6 

-3 

21 

ll 

lbub.bu 

319 

_ 

269 

35b 

256 

Comparison  of  calculated  dipole  matrix  elements  with  experimental  and  rigid 
rotor  values:  (ulo)  band,  147o  cm*^  (mostly  P  branch). 

■  experimental  value*^' 

S..,,.  *  present  results 

Cd  l  c 

( lb ) 

SD  ■  previous  results 

O 

S  ■  rigid  rotor  value 

rigid 

j  All  values  in  (debyes)2 
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TABLE  5  (concluded) 


7-2 

! 

1^54 .99 

I 

lbl 

114 

- - 

133 

466 

J-1 

4j 

l^bb.27 

.16 

.13  | 

.20 

1.06 

7 -4 

B-2 

1961.19 

bl.U 

l 

37.3 

66.9 

1 

161. 

*2 

7  2 

1900.27 

236 

140 

216 

72U 

4-2 

b2 

197b. 19 

.96 

.66 

1.06 

5.36 

i _ !l_ 

i _ 

190B.41 

220 

14u 

212 

7  00 

Comparison  ot  calculated  dipole  matrix  elements  with  experimental  and  rigid 
rotor  values:  (olu)  Band,  1900  cm  ^  (mostly  R  branch). 

S  *  experimental  value'^* 

exp  r 

S. . i .  *  present  calculated  value 

Cd  I  c 

(  lb  ) 

S  *  previous  results 

b 

S  4  ^  *  matrix  element  for  rigid  rotor 

rigid 

All  values  in  (debyesr 
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OR  WEAK  LIN 


Comparison  or  calculated  dipole  moments  with  calculated  and  rigid  rotor 
values:  ( olo )  oano,  <!35u  cm-1  (R  branch).  (No  experimental  results  avail 


S  ■  present  cal cul ate<Lya1 ue 

Sy,c  «  previous  results'18’ 

.  ■  rigid  rotor  value 

All  Values  in  (debyesr 
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IIIC.  Line  Broadening  and  Shifts 

We  assume  that  broadening  is  dominated  by  collision  processes.  In 
particular,  we  calculate  the  broadening  and  shifts  of  H,0  due  to  collisions 
with  N0  and  On,  the  major  atmospheric  components.  Anderson  theory  discussed 
above  is  used.  All  that  remains  is  to  choose  the  form  of  the  interaction 
between  the  molecules  which  mediates  the  broadening  process. 

In  the  case  of  ^  collisional  deactivation  by  we  found  that  long 
ranged  electromagnet ic  interactions  were  very  ineffective  deactivation  agents. 
The  long  ranged  interactions  are  highly  resonant,  leading  to  little  energy 
transfer  for  the  case  aE  v>  760cm’1.  Only  through  the  mechanism  of  hard  core 
repulsive  forces  can  significant  VV  exchange  occur. 

To  determine  the  amount  of  line  broadening  in  Anderson's  theory,  the 
transition  probabilities  to  all  possible  neighboring  states  are  sunned. 
Consequently  the  nearest  neighbor  transitions  will  dominate.  To  deal  with 
pure  rotational  transitions  with  relatively  small  energy  exchange,  long  ranged 
interactions  are  sufficient.  In  general,  hard  core  collisions  lead  to  large 
energy,  vibrational  exchanges,  while  long  range  interactions  cause  small 
energy  rotational  exchanges.  The  collis'onal  broadening  and  shifts  of 
spectral  lines  are  dominated  by  the  long  ranged  electromagnetic  interaction. 
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The  lowest  order  electromagnet ic  interaction  between  H^O  and  N?  or  02 
is  d i pol e-quadrupol e.  For  this  case,  it  is  well  known  that  the  response  func¬ 
tion  r(w)  resulting  from  the  integral 


r(w)  -  |£  ^  e',wt  V(t)  dt|2 


(3-9) 


is  a  sun  of  products  of  Kelvin  functions. 

Consequently,  the  collisional  line  broadening  can  be  readily  deter¬ 
mined  using  (2-14)  and  (2-16).  For  dipole-quadrupole  interaction,  (2-14)  for 
the  level  |vJr>  of  H90,  J0  of  N,  becomes 

c  c  c 


I <  v J  r | u | v'J“V>|‘  | <J . I Q| Jp> | 2  C2(J.2J,;00) 


•  h^v  btt  ^  '  -i-i’  -  •  i 


v  J  t  ,J, 


where 


(3-10) 


G(wt)  «  ^  [  (3+u) !  (3-u)!]'1  (wt)6  kjlwtl) 


(3-11) 


101 J2>  *  quadrupole  moment  for  or  0^ 

C(J^2J2;00)  ■  Clebsh  Gordon  coefficient 

W  -  energy  defect  «  E-,.  E,  -  E  <  -  E  , 

v  J  r  “J' 
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For  the  case  of  vibrational  deactivation  we  neglected  the  rotation 
entirely,  because  the  vibrational  exchange  was  orders  of  magnitude  greater 
than  the  rotational  energy  change.  For  line  broadening  all  the  possible 
transitions  are  sunned  over,  and  the  smallest  energy  exchanges  dominate.  The 
smallest  energy  exchanges  occur  in  purely  rotational  transitions,  and  thus  the 
rotational  energies  can  no  longer  be  neglected.  Since  quadrupole  selection 
rules  are  obeyed,  the  sum  in  (3-10)  over  J,  ranges  over  a  maximum  of  five 
values  for  each  J^,  so  this  additional  complication  is  not  serious.  When  the 
initial  states  J,  of  the  broadener  are  averaged  over  as  in  (2-16),  the  sum  is 

_9 

truncated  for  those  J,  Such  that  the  thermal  population  is  less  than  10  of 
the  ground  state.  This  occurs  for  J,  <  100  for  most  temperatures . 

To  calculate  the  response  function  (2-15)  for  the  collisional  line 
shifts,  the  Hilbert  transform  of  (3-11)  must  be  evaluated.  Instead  of  working 
directly  with  (3-11)  for  this  purpose,  we  used  the  least  squares  fit 


G(x)  *  e‘2x  (.1777  ♦  .35l9x  ♦  .3877x2  ♦  .2759x3  ♦  .1764x4  ♦  .1395x5) 

(3-12) 


Denote  Ik  the  Hilbert 


transform  of  |w|k  e'^w^T; 


(3-13) 


Then  1^  can  be  expressed  in  terms  of  exponential  integrals: 


2w  r 


-2w„  r 


^  [+e  0  E1(2w0r)+e  0  Ei (2wqt)3 


W„T  2w„I 


-2w0t 


+  ~k  [”e  °  E1(2wot)  +  e  6,(2-003 


w  r 

l2  "  '  77" 


(w  tj  2w  t  -2w  t 

*  Oe  0  E1(2wqt)  ♦  e  0  E^i)] 


woT 


:w  t 


w  r )  2w  r  _ 

♦  — rr  C-e  0  E1(2w0t)  ♦  e  '~0'  Ei(2w0t)] 


Wn<  (Wn03  !«  t)4  _  2wrt 


-2w, 


l4’-zr  -TT-  °  El<2V>+e  °  E1<2V» 


3w  T  (w  r)^  (w  t)^  2w  t  -2w  t 

“37 - TT  *  ~~7T  [*e  E,  (2w00  ♦  e  °  £,(2^03 


In  the  limit  of  large  «0 t. 


lk  '  TT 


1  V  (k  ♦  2s)! 

2k  w„t  Zj  (2w„  02* 


(3-15) 


S-0 


Thus  the  line  shift  response  function  drops  off  only  as  1/w,  while  the  line 
broadening  response  function  falls  off  exponentially.  However,  the  line  shift 
response  function  is  an  odd  function  of  wQ,  leading  to  a  cancellation  between 
terms  of  positive  and  negative  energy  defects.  Therefore,  even  though  *(w) 
has  a  much  wider  range  than  r(w),  the  level  shifts  and  widths  resulting  from 
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(2-16)  are  roughly  the  same  size  since  r(w)  is  an  even  function  and  »(w)  is  an 
odd  function  of  w.  The  velocity  averaged  >  and  r  as  functions  of  x  *  wb/M/2KT 
are  pictured  in  Figure  2.  Note  that  though  r  is  larger  for  small  w,  it  de¬ 
creases  exponentially  while  •  is  decreasing  only  as  1/w. 

Finally,  we  must  average  over  v  and  b.  The  integral  over  v  was  done 
mnerical  ly.  Since  r(b)  *  •  for  b  •  0,  the  assumption  that  r  is  small  is 
invalidated  in  that  region.  We  used  an  averaging  process  due  to  Sharma  and 
Srau'1^  to  eliminate  the  necessity  of  numerically  integrating  over  b  using 
the  full  exponential  form  for  (2-16).  We  assune  that  r(b,v)  is  given  by  the 
fol lowing: 


r(b)  *  r(b)  r  l/b6  b>d 


r(b) 


r(o) 


■i 


C r ( d)  -  r(o) ] 


b<d 


d  *  Lennard-Jones  a 


(3-16) 


To  obtain  the  relation  for  r(o),  we  recalculate  (3-9)  for  head-on  collisions 
only.  A  similar  result  is  obtained,  with  (3-12)  replaced  by 


G0(‘) 


J?(«) 


-x/2 


(.1100  ♦  . 0392x  -  •0295x‘"  ♦  .0064x' 


-  ,0006x<:) 

(3-17) 
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Dipole-Quadrupole  Response  Functions 
for 

Widths  and  Shifts  of  Levels 


’’he  Hilbert  transform  can  be  evaluated  as  before.  Averaging  (3-16)  over  b,  we 
obtai  n 


r(b)  bdb  -  |  r(o,v)  +  £  r(d.v)  (3-17) 

The  averaging  integral  over  v  was  done  numerically,  leaving  a  line 
broadening  response  function  and  a  line  shift  response  function  depending  only 
on  energy  defect  w.  Finally,  the  averaging  over  the  N0  initial  state  was  done 

C 

assuning  a  Boltsnan  distribution: 


P( Jj)  *  ( 2 J j ■*"  1 )  e  3J,*vJl  +  1)  /  £  (2J*1)  e*BJ(J+1)  (3-13) 

J 


For  H^O,  we  used  the  dipole  moment  operator  given  by  3.6.  The  quadru- 
pole  moments  and  rotational  constants  for  N.,  and  0~  are  given  below. 


N^,  and  0^,  Molecular  Parameters: 


Q 

b 

( x  esu) 

(cm-1 ) 

H2 

1.7 

2.U1 

°2 

l.U 

1.45 

In  Tables  a  ano  9  we 

give  the  resulting 

level  hal f  widths 

ano  shifts 

due  to 

and  0^  at  T*juu0. 

Note  that  to  obtain 

the  line  halt  width,  the  two 

1  evel 

naif  widths  are  added 

.  Conversely,  to 

obtain  the  line 

shift,  we 

subtract  tne  lower  level  line  shift  fran  the  upper  level  line  shift.  Conse¬ 
quently,  the  shirts  are  generally  an  order  of  magnitude  less  than  the  widths. 
In  fact,  the  shifts  are  so  small  for  atmospheric  pressures  (»  .ul  cm"^),  that 
tney  are  less  than  the  accuracy  of  our  model.  In  Taole  2u  we  compare  the 
halt  widths  with  other  theories  ano  experiments.  The  results  are  in  general 
very  good,  with  differences  attributable  mainly  to  the  approximate  impact 
parameter  averaging  process  used.  To  see  this,  we  need  only  note  that  the 

_4 

widths  depend  sensitively  on  the  cutoff  impact  parameter  as  d  •  A  numerical 
impact  parameter  averaging  may  De  used  with  the  exponential  form  of  the 
propaoility,  but  this  leads  to  a  factor  of  ten  increase  in  computation  time. 
If  additional  accuracy  is  needed,  this  modification  can  easily  be  done,  and 
the  resulting  calculation  will  consume  approximately  1  minute  of  CDC  7buu  time. 
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2.52 

1.2a 

2.5tt 


1.97 


TAbLE  9.  CALCULATED  WIDTHS  AND  SHIFTS  OF  ROTATIONAL  LEVELS  OF 
H20  DUE  TO  COLLISIONS  WITH  02>  IN  MHz/T 


TABLE  lu.  COLLISION  BROADENED  LINE  WIDTHS 


H20  BROADENED  BY  02  AND  N^  COMPARISON  WITH  OTHER  THEORIES  AND  EXPERIMENT 


1  ! 

Broadened  (mhz/T)  J 

Broadened  ( MHz/T ) 

Line 

Theory 

Present 

T 

Experiment  ( 

jlheory * 

Present 

Experiment 

lul"llu  j 

2.67 

1.14 

( 

I 

4.14 

1.68 

*o2”zll 

2.6b 

1.78 

j 

4.09 

2.47 

^2u“J13 

2.06 

1.47 

2.7 

4.  lu 

2.24 

3.8 

J21"4I4 

2.6* 

1. 62 

4.16 

2.64 

2.6* 

l.U* 

4. 14 

1.67 

422_bl6 

2.73 

2.06 

4.20 

3.2* 

2.74 

.7* 

4.23 

1.41 

441'532 

2.63 

l.lu 

4.06 

1.71 

b2J*°lo 

2.74 

2.20 

2.6 

4.23 

3.62 

4.1 

2.6* 

.64 

4.14 

1.03 

|  35i'°42  | 

2. 7u 

.66 

4. 16 

1.16 

' 

1 1 IU .  oand  Model 


The  line-by-line  band  model  is  based  on  a  list  ot  line  frequencies, 
strengths,  ano  widths.  First,  a  table  of  collisional  widths  and  shifts  for 
the  specific  temperature  and  composition  is  formed.  Then,  all  the  possible  VR 
transitions  are  determined,  the  intensities  ot  each  line  calculated,  the 
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widths  determined  by  adding  the  widths  of  the  upper  and  lower  level,  and  the 
shifts  by  subtracting  the  two  level  shifts. 

In  Appendix  1  we  give  a  list  of  the  water  vapor  lines  and  strengths 
for  the  mF  and  DF  regions  of  the  spectrum  2ttuu  to  4bUU  cm'*.  The  half  widths 
ana  line  shifts  for  each  line  for  a  specific  concentration  of  and  0^  can  be 
determined  from  tables  7  ana  y.  Since  these  tables  are  determined  tor  the 
ground  state,  the  calculation  should  be  repeated  for  each  band,  since  the 
a symmetry  varies  with  the  degree  of  vibrational  excitation.  However,  to  a 
first  order  of  approximation  the  results  for  the  ground  state  can  be  used  for 
any  band. 


The  calculation  of  the  wiatns  and  shifts  for  each  band  can  easily  be 
adopted  into  the  program  which  calculates  the  lines  ana  strengths  to  give  all 
the  parameters  of  a  band  model  for  each  temperature,  pressure,  and  composi¬ 
tion.  The  exact  f"rm  of  such  a  band  modeling  program  obviously  would  depend 
on  its  purpose. 

C.  Description  of  Computer  Codes 

I .  I ntroductlon 

The  present  codes  consist  of  a  set  of  12  modules  whose  interrelations 
are  Indicated  in  Figures  3  and  4  tor  rate  calculations  and  band  modeling.  In 
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Fioure  3 

Code  structure  for  rate  calculations 
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general,  data  is  transferred  from  module  to  module  by  storing  it  on  disc.  A 
brief  description  of  each  module  follows: 

1 1 .  Description  of  Codes 

(a)  Model 

This  program  has  been  documented  elsewhere.  Its  purpose  is  to 

calculate  the  non-diagonal  rotation  matrices  A,  B,  and  C.  It  does  so  by  vary¬ 
ing  the  vibrational  potentials  to  fit  observed  band  origons.  Output  quanti¬ 
ties  are  the  band  origins  and  the  rotation  matrices.  If  needed,  the  Morse 

oscillator  matrix  elements,  the  vibrational  wave  functions,  and  the  rotation 
matrices  are  stored  on  disc  for  later  use.  The  vibrational  wave  functions  are 
the  coefficients  of  an  expansion  in  a  set  of  Morse  oscillator  functions. 

(b)  STAND 

STANO  consists  of  two  utility  programs,  HDIAG  and  STEP  IT.  HD1AG 
diagonalizes  syrnetric  matrices  and  calculates  the  necessary  transformation  ma¬ 
trices.  STEPIT  is  the  fitting  routine.  The  parameters  to  be  varied  are 

stored  in  the  array  Z([),  and  the  sum  of  the  squares  of  difference s  to  be 

minimized  is  calculated  by  the  main  program  routine  SHSQER .  MOOEL,  BAND,  and 

LINE  must  be  loaded  with  the  utility  programs  in  STAND. 
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(C) 


ROTATE 


This  program  calculates  the  roto- vibrational  energy  levels  using  the 
rotational  matrices  stored  on  disc  by  MODEL.  If  needed,  the  energy  level  as 
well  as  the  roto- vibrational  wave  functions  in  the  form  of  coefficients  of  an 
expansion  in  vibrational  wave  functions  and  symmetric  top  wave  functions  are 
stored  on  disc.  These  coefficients  take  up  large  amounts  of  storage  space, 
but  storage  is  convenient  since  this  data  is  used  by  three  other  programs  for 
different  purposes. 

The  wave  functions  resulting  from  MODEL  are  used  in  WANG  where  the 
matrix  elements  of  the  rotational  piece  of  the  Hamiltonian  in  the  basis  of 
Wang  states  are  calculated.  The  resulting  matrix  is  diagonalized  to  yield  the 
roto- v ibrat i onal  energies  and  transformation  matrices.  The  C2y  symmetry  is 
taken  fully  into  account.  The  roto-vibrational  wave  functions  are  then  repre¬ 
sented  by  the  columns  of  the  transformation  matrices  U,  and  are  just  the  co¬ 
efficients  of  an  expansion  in  a  basis  space  consisting  of  a  product  of  vibra¬ 
tional  wave  functions  and  Wang  symmetric  top  wave  functions. 

(C)  IFN 

This  program  calculates  the  fourier  transform  response  function  of 
any  multipole-multipole  electrtjnagnetic  interaction  in  the  classical  path 
approximation.  As  is  well  known,  the  result  is  given  in  terms  of  a  sum  of 
products  of  Kelvin  functions.  In  addition,  the  response  function  for  a  head 
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on  collision,  with  impact  parameter  0*0,  is  calculated.  The  results  can  be 
used  for  V,  VR,  or  RR  rates. 


(e)  QFN 

This  routine  takes  the  ordered  energy  levels  placed  on  disc  by  ROTATE 
and  evaluates  the  statistical  Q  function  given  below  as  a  function  of  tempera¬ 
ture: 

Q  *  gjir(2J*l)e"Ejr/KT  (2-1) 

(f)  RTDAT 

This  program  evaluates  the  matrix  elements  of  an  arbitrary  potential 
V(R,»,»)  between  VR  states  whose  energy  difference  lies  in  a  particular  spec¬ 
tral  region  and  places  the  results  on  disc.  The  potential  is  expressed  as  a 
series  of  spherical  harmonics,  and  the  result  is  found  using  the  known  matrix 
elements  of  speherical  harmonics  between  symmetric  top  wave  functions.  The 
full  roto-v ibrat • onal  wave  functions  found  by  ^ODEl  and  ROTATE  are  used. 

(g)  RATE 

RATE  calculates  transition  probabi 1  it ies  induced  by  an  arbitrary 
potential  V  by  confining  the  matrix  elements  of  V  found  by  RTDAT  with  an  appro¬ 
priate  response  function.  The  result  is  given  as  a  function  of  temperature. 
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BAND 


(h) 

To  determine  the  dipole  matrix  element  given  by  (3.1),  this  program 
varies  the  in  order  to  fit  experimental  integrated  band  intensities.  The 

vibrational  wave  functions  resulting  from  MODEL  are  used,  and  the  fitting 
routine  is  STE PIT. 

(i)  LINE 


LINE,  TABLE,  and  WIDTH  all  use  the  roto- vibrational  wave  functions 
found  by  ROTATE  to  calculate  dipole  moment  matrix  elements.  The  dipole  moment 
operator  is  gven  by  the  relation 

“  ■  £  “Ijk  ql  4  q3  <2-2> 

where  u,,^  have  the  usual  dependence  on  rotational  coordinates  and  are 

vibrational  coordinates.  LINE  calculates  the  dipole  transition  strength  for 
up  to  500  H2O  absorption  lines,  which  are  given  as  input.  An  option  is  to  use 
the  fitting  routine  STE P I T  to  vary  the  parameters  in  order  to  acheive  a 

fitting  of  experimental  line  strengths. 


(j) 


WIDTH 


This  program  uses  the  Anderson  theory  of  collision  broadening  to 
calculate  the  widths  and  shifts  of  H^O  levels  caused  by  collisions  with  homo- 
nuclear  molecules  such  as  N^  and  0V  The  dipole  quadrupole  moment  long  ranged 


6T 


interaction  is  presumed  to  dominate  the  process.  The  probabilities  of  virtual 
transitions  to  all  possible  neighboring  levels  are  calculated  and  summed  up  to 
yield  both  the  level  width  and  shift.  Each  individual  probability  is  assimed 
to  be  small  enough  so  that  the  non-exponential,  single  quantum  exchange  forn 
of  the  transition  probability  can  be  used.  In  this  case,  an  approximate  aver¬ 
age  over  the  impact  parameter  b  can  be  used,  wh.ich  eliminates  the  need  for  a 
numerical  integration.  This  approximate  averaging  was  also  used  in  our  colli- 
sional  rate  calculations  and  gives  only  20-40i  accuracy,  but  saves  quite  a  bit 
of  computer  tine.  If  accuracy  is  needed,  the  exponential,  multi-quantum  ex¬ 
change  form  of  the  transition  probab i 1 i t ies  together  with  numerical  integra¬ 
tion  over  impact  probabilities  should  be  used.  The  line  widths  and  shifts  in 
units  of  cm" V  atm-*  are  given  for  each  level.  Those  for  the  highest  J  states 
included  are  wrong,  because  the  virtual  transitions  to  J+l  were  not  included. 


(k)  I  FUNS 

In  this 
sponse  functions 
transform  of  tne 

G(x)  *  e 

where  x  »  wb/ v 

and  G ( x)  •  | 


specialized  program  the 
for  line  shifts  and  line 
interaction  notential  is 


i 

K 


j (t)dt| 


2 


velocity  and  impact  parameter  re¬ 
widths  are  calculated.  The  fourier 
fit  to  the  form 

(2-3) 

(2-4) 

(2-5) 
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Then  the  line  shift  response  function  is  given  in  terns  of  exponential  inte¬ 
grals.  In  addition,  the  head  on  response  function  is  determined.  Both  are 
averaged  over  velocity,  and  the  approximate  impact  parameter  averaging  de¬ 
scribed  in  part  8  is  performed.  The  results  are  used  in  WIDTH,  where  the 
response  functions  are  determined  by  interpolating  between  values  calculated 
by  I FUNS. 

(1)  TABLE 

TABLE  uses  the  energy  levels  from  ROTATE  to  create  a  spectrum  of  all 
possible  absorotion  lines.  The  dipole  strengths  are  calculated  assuming  a 
dipole  operator  given  by  3.1,  and  the  results  are  given  in  the  form  of  a 
table,  ordered  by  1'ne  frequency.  Both  the  temperature  dependent  line 
strength  and  the  dipole  matrix  element  is  printed. 
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APPENDIX  1 


Table  of  Absorption  ' ines  from  28UU  cm"*  to  4800  cm 


Explanation  of  Table: 


Col  jnn: 


1.  lower  state  vibrational  quantum  nunoers  (v^  v.,  v^) 

2.  upper  state  vibrational  quantum  numbers  (v,  v3) 

3.  J  of  lower  state 

4.  t  of  lower  state 

5.  J  of  upper  state 

c.  t  of  upper  state 

7.  Frequency  of  line  In  cm’ 

8.  Energy  of  lower  level  in  cm" 

9.  tine  strengtn  in  cm’Vgm  cm* 

10.  Transition  dipole  moment  in  (debye)  xlo  x  (2J*1) 

T  «  300°K 

bands  included  ■  (OOu),  (UlO).  (100),  (110).  (020).  (120),  (OOl),  with  J  <  10 
Centrifugal  distortion  included  for  all  but  the  (001)  band. 

Coriolis  rnteraction  neglected. 

Oipole  moment  operator: 

y  ■  1.8519  ♦  .4580  qj  -  .  Toll  q^  -  .4*74  q^ 

-  .7890  q, q^  ♦  .0143  q*  ♦  1.4572  q3 


The  line  strength  at  any  temperature  can  be  obtained  from  Column  10  using  the 
following  relation: 
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